The origin of magnetic freezing in the pyrochlore Y2M02O7 
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We investigated the nature of the spin glass-like phase transition in the geometrically frustrated 
pyrochlore lattices Y2M02O7 using the local probes nuclear and muon magnetic resonances, and the 
field-dependent long range probes x-ray and neutron scattering. The long range probes indicated 
that Y2M02O7 does not undergo any global symmetry changes, even in a field of 6 T. In contrast, 
the local signal indicates a lattice distortion close to the critical temperature. The nuclei show 
at least two inequivalent Y sites, and the muons show sub-linear line broadening as a function of 
moment size, over a wide temperature range. The conclusion from all the measurements is that 
even in high field, the distortion of Y2M02O7 takes place within the unit-cell, while its global cubic 
symmetry is preserved. Moreover, the muon result clearly indicates the presence of magneto-elastic 
coupling. 

PACS numbers; 75.50.Lk, 75.10.Nr 



The Heisenberg model on the geometrically frus- 
trated pyrochlore lattice has a macroscopically degen- 
erate ground state, and the standard degeneracy-lifting 
mechanism of thermal or quantum fluctuation does not 
seem to remove it. Yet, with only one exception, the 
family of compounds based on Tb2Ti2 07 \^\^, all py- 
rochlores freeze (at least partially), that is, one state out 
of many is selected. In some cases, like the spin ice, 
this is due to long-range interactions [5] and single ion 
anisotropy In others, the freezing occurs even with- 

out anisotropy. In these cases magneto-elastic coupling 
might be responsible for the degeneracy lifting; the lat- 
tice distorts to relieve the frustrated interactions. Such a 
distortion might lead to a cubic-to-lower-symmetry struc- 
tural transition 8]. This kind of frustration-driven dis- 
tortion has been previously suggested as the main freez- 
ing mechanism for several ^2^207 pyrochlores and 
Cr spinels 11 , 13] , and was considered theoretically 13 - 
Il6j . However, lattice distortions and symmetry changes 
as a function of temperature are a common feature in 
solids, even without magnetic interactions. Therefore, it 
is not yet clear whether: (I) the magnetic interactions 
drive the distortion; or (II) the distortion takes place 
because of electrostatic interactions, and the magnetic 
properties, such as freezing, follow. Clarification of this 
point is crucial for the understanding of the spin Hamil- 
tonian, and therefore the ground state and excitations in 
pyrochlores. Field-dependence experiments can provide, 
in principle, answers to these questions. 

One case is the pyrochlore Y2M02O7 (YMO), which 
crystallize into a cubic structure with Fd3m symmetry 
17| . Magnetically, it has a Curie- Weiss temperature of 
—200 K and freezes with spin-glass characteristics [is'l 



at Tf — 22.5 K 'lo]. In particular, magnetization mea- 
surements indicate a large difference between zero-field- 



cooled and field-cooled magnetizations [20|, |21|. This 
glassiness is unexpected if the chemical structure is per- 
fect because it is believed that a spin-glass state emerges 
when frustration and disorder coexist. A possible so- 
lution to this dichotomy came from more detailed mea- 
surements. X-ray absorption fine structure spectroscopy 
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FIG. 1. (Color online) The raw muon depolarization data 
taken at T = 24K, at two transverse fields. The error bars 
represent statistical errors. The inset shows the bulk suscep- 
tibility X = M/ H versus temperature in zero-field cooled and 
field cooled conditions taken at a IkG field. 
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FIG. 2. (Color online) The TF-^iSR relaxation rate, I/T2, 
versus the magnetization, M at T = 24K. The field jxoH is 
an implicit parameter. Hollow circles are taken from Ref. [l^] • 
The red (solid) line is a linear fit to low Ad data. The black 
(dotted) line is a fit to A • M'' where /3 = 0.66(12). The inset 
shows l/Ta* versus M at T = 100 K. 



(XAFS) in zero field show evidence of positional disorder 
of the Mo ion In contrast, neutron pair distribution 
function measurements assigned the distortion to the 01- 
Y bond [2^. Local magnetic probes, which are coupled 
to the spin system and operate in a field, such as ^^Y 
nuclear magnetic resonance (NMR) [9f and muon spin 
relaxation (/iSR) [lo| . indicate a distortion of the spin- 
bearing ion, again the Mo. Despite the controversy on the 
distorted ion, all experiments suggest the lattice is tak- 
ing part in the magnetic freezing. However, the presence 
of an applied field in the resonance measurements and its 
absence in scattering measurements make the comparison 
difficult. Moreover, neither experiment can shed light on 
the magneto-elastic coupling issue discussed above. The 
motivation of this work is to fill in the gap and perform 
both local resonance and long range scattering measure- 
ments under the same conditions, and to use the field as 
a probe of magneto-elastic coupling. 

We performed five different experiments on YMO, 
which were carried out well above and close to Tf. (i) 
High transverse field (TF) and longitudinal field (LF) 
/iSR, which extends previously low-field data Eol. HJ]- 
(ii) ^^Y NMR where we extend previous data [9| to the 
helium range, (iii) Field-dependent high resolution x- 
ray powder diffraction, (iv) Field-dependent neutrons 
diffraction, which extends the previous zero field (ZF) 
measurement 



171 . |23|. (v) Bulk magnetization using a 



superconducting quantum interference device (SQUID). 

Polycrystalline samples of YMO were prepared accord- 
ing to Ref. 19]. The inset of Fig. [T] shows the zero-field- 
cooled (ZFC) and field-cooled (FC) magnetization mea- 



FIG. 3. (Color online) The exponent /3 versus the Tempera- 
ture. P is extracted from I/T2 versus the magnetization M 
fits to I/T2 = AM^ . Line is a guide to the eye. 



surements. The ZFC curve shows a distinctive maximum 
indicating the spin-glass transition. Tf and Curie- Weiss 
temperatures extracted from these measurements are in 
agreement with previous reports. In addition magnetiza- 
tion vs. applied field up to 5 T at various temperatures 
was also recorded and will be used below. 

^SR measurements were performed at TRIUMF, 
Canada on the HiTime and Helios spectrometers in the 
M15 surface muon channel. Figure [1] depicts the raw 
muon polarization taken at a constant temperature, T = 
24 K, at two transverse fields /ioi?TF = 0.4 and 5 T. 
The data is shown in a rotating-reference frame [2^ of 
Mo^^RRF = Mo^^TF ~ 0.02 T. The TF relaxation rate in- 
creases with increasing fields. The LF (not shown) is 
at least an order of magnitude smaller than the TF re- 
laxation rate at any field and temperature. This means 
the muon depolarization is mainly due to static field in- 
homogeneities and the contribution of the out-of-plane 
depolarization can be neglected. We found that the /xSR 
TF asymmetry is best described by 



PtfW = Poexp 



-^t/T*^cos{ujt + ^) (1) 



where Pq is the initial polarization and w = 7^iJTF- The 
fits are represented by the solid line in Fig. [TJ 

In Fig. [2] we plot the transverse field relaxation rate 
I/T2* at a given field versus the magnetization M mea- 
sured at the same field, at temperature fixed at 24 K, 
slightly above the spin-glass freezing. We also show the 
field, iJ, as an explicit parameter on the upper abscissa. 
As the field increases the relaxation rate I/T2* increases. 
However above M ^ IQ-^ Am^, is no longer a 

linear function of M . This is not the case at higher tem- 
peratures. For comparison, in the inset of Fig. [5] we plot 
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FIG. 4. (Color online) The NMR spectra at 25 < T < 300 K. 
Inset: the NMR spectrum at T = 25 K. The line is a fit to 
a powder spectrum with three sites. The arrows point to the 
main peak in the powder spectrum of each site. 



I/T2* vs M taken at 100 K, again with H as an implicit 
parameter. At this temperature the relaxation rate is a 
linear function of M . 

In a system with quenched disorder it is expected that 
I/T2* oc M 0, m as is indeed observed at 100 K. How- 
ever, when T approaches T/, the data deviates from this 
linearity and 1 /Tj* grows more slowly than M as depicted 
in Fig. 13 Fits to l/T^ = AMl^ are also shown in the fig- 
ure, giving (3 = 0.66(12) and /3 = 1.00(2) for T = 24 
and 100 K respectively. It is found that the transition 
from /3 = 1 to very low values are spread over a wide 
temperature range (Fig. [3]). 

One possibility is a site-dependent spin polarization 
due to impurities [27| . However, this can be ruled out 
by the undetectable amount of disorder in the sample. 
A more likely scenario is that the disorder in YMO is 
not quenched and at low T it varies with the field. In 
fact, the sub- linear behavior suggests that as the field 
(and M) increases the lattice is more ordered. This is 
the main finding of this work, which clearly points to the 
presence of a magneto-elastic coupling. 

More evidence for temperature dependent disorder 
comes from NMR which extends previous measure- 
ments down to 25 K. Such low-T measurements were 
made possible by using a high pressure cell where higher 
RF power can be delivered to the sample. At each 
temperature we obtain the complete NMR spectrum by 
sweeping the external field, H^xti at a constant applied 
RF frequency / = 16.44 MHz. In each field we used 
the spin-echo sequence (7r/2 — tt pulses) and recorded the 
echo signal. In Fig. 2] we present the spectra taken at 
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FIG. 5. (Color onhne) The NMR shift, Ky, versus the bulk 
susceptiblity x a-nd the temperature as an implicit parameter. 
In the inset, a theoretical powder spectrum of a nuclear spin 
1/2 is depicted. 



temperatures between 300 K and 25 K. The width of the 
300 K spectrum extends over 0.1 T, whereas the width 
of the 25 K spectrum extends over 1 T. This broaden- 
ing results in low intensities at each applied field upon 
cooling. Due to this broad line at low T we gave up on 
high resolution NMR, as in Ref. ^ , and concentrated on 
the gross features of the spectrum. The most noticeable 
feature in the T = 25 K spectrum is the clear appearance 
of two peaks, with a hint of a third one. This suggests 
that out of the many different Y sites existing at high 
r[9], only few are being picked as T is lowered. 

To understand this NMR spectrum we look at the 
spin 1/2 ^^Y nuclear spin Hamiltonian, which can be de- 
scribed as. 



^hjl • (1 + K) • He 



(2) 



where K is the NMR shift tensor, I is the nuclear spin op- 
erator, and Hext is the variable applied field. In powders, 
the principal axes of K are randomly oriented relative to 
Hext- Therefore, the magnetic resonance spectrum is an 
average over all possible orientations. A theoretical pow- 
der averaged NMR line [28[, for a single site, is depicted 
in the inset of Fig. [5] where Ha = 27r//[7(l -I- Ka)] and 
a — x,y,z for the three directions. This theoretical spec- 
trum demonstrates that a single site, with a single set of 
Kx, Ky, Kz, could give rise to only one peak even under 
powder averaging. The existence of two (perhaps three) 
peaks in the spectrum is a result of a lattice deformation 
leading to inequivalent Y sites. 

In order to study the temperature dependence of 
the shift, we use the powder spectrum convoluted with 
Lorentzians to fit the NMR spectra. Such a fit is demon- 
strated by the solid line in the inset of Fig. S) At high 
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temperatures (T > 250 K) the shift is very smah; at in- 
termediate temperatures, 50 < T < 200 K, two different 
sites were needed to fit the data (see Fig. and finally, 
at low enough temperatures, T < 50 K, three sites were 
assumed. In Fig. [5] we plot the shift K^, i for each site, 
versus the bulk susceptibility, which was extrapolated to 
fioH = 7.8 T from the FC magnetization measurements. 
Temperature in this figure is an implicit parameter. As 
X increases the shift for each site also increases. How- 
ever, the dependence between shift and susceptibility is 
not linear. When the disorder is quenched one expects 
K (X X- This proportionality is violated close to T/, in- 
dicating that the lattice degrees of freedom are active as 
T ^ Tf. It should be pointed out that unlike in /iSR, 
in NMR it is impossible to vary the field over a wide 
range and NMR cannot be used to address the question 
of magneto-elastic coupling. In contrast, the presence of 
multiple sites well above Tf suggest that spin correlations 
are sufficient to distort the lattice. 

We also searched for field effects in high resolution X- 
ray scattering. The X-ray powder diffraction experiments 
were conducted at the APS Argonne National Laboratory 
on the 11-ID-C beamline. A high-energy, 115 keV x-ray 
beam with a high-resolution analyzer was used with a 
6 T magnet. In order to dismiss any grain orientation 
with the field, GE varnish was applied. In Fig. [DJa) we 
plot the most intense cubic (440) peak with and without 
an applied field at 60 K and 25 K. The peaks are field- 
independent in both shape and intensity. Needless to 
say, no peak splitting or new peaks were found when the 
field was applied and moreover the peaks are resolution 
limited. This rules out any global structural transition 
due to temperature or field. 

Finally, similar experiments were performed using neu- 
trons which are more sensitive to scattering from oxygen, 
but have lower resolution. These experiments were done 
on the BTl powder diffractometer at NIST, Gaithers- 
burg, USA, with a field up to 6 T perpendicular to the 
scattering plane. Data were collected at the same tem- 
peratures, with a neutron energy of E = 34.5 meV. In 
Fig. E^b) we plot the (440) Bragg peak. As with the x-ray 
picture, no apparent difference is revealed between the 
measurements with and without the field. These peaks 
are resolution limited as well. However, the small dif- 
ference in intensity at T = 25 K between the 6 T and 
zero field measurements is real and reproducible. It is 
probably due to very small structure-factor or magnetic 
form-factor changes induced by the field. 

To conclude, both local probes unambiguously indicate 
a lattice deformation takes place as T approaches Tj. In 
addition, the magnetization dependence of the muon spin 
relaxation rate shows that lattice deformation is affected 
by magnetic field, therefore pointing to magneto-elastic 
coupling. Long range scattering measurements fail to 
detect global changes in the lattice parameters upon ap- 
plication of the field. Since the distortion is found locally 
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FIG. 6. (Color online) Field and temperature dependence of 
the (440) Bragg peak from (a) X-ray and (b) neutron scat- 
tering. The error bars representing statistical errors of ±la 
are smaller than the symbols. The peaks in (a) and (b) are 
resolution limited. 



by resonances, but not globally by scattering, it must be 
within the unit cell. It affects, at most, the structure- 
factors on the one hand, and the hypcrfinc coupling be- 
tween Y nuclei and Mo spin on the other. Our work indi- 
cates that magneto-elastic coupling is part of the freezing 
process of YMO and provides a simple way to detect it 
in other magnets. 
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